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We prove that taking correctly into account the lepton flavour dependence of the CP asymmetries 
and washout processes, it is possible to obtain successful thermal leptogenesis from the decays of the 
second right-handed neutrino. The asymmetries in the muon and tau-flavour channels are then not 
erased by the inverse decays of the lightest right-handed neutrino, Ni. In this way, we reopen the 
possibility of “thermal leptogenesis” in models with a strong hierarchy in the right-handed Majorana 
masses that is typically the case in models with up-quark neutrino Yukawa unification. 


The measure of nonvanishing neutrino masses and mixings and the use of the seesaw mechanism to explain the 
smallness of these masses has made leptogenesis [lj one of the most promising mechanisms to generate the observed 
value of the baryon asymmetry of the universe. The seesaw mechanism requires the presence of heavy right-handed 
neutrinos with complex Yukawa couplings that can generate a lepton asymmetry through their out-of-equilibrium 
decays. This “lepton” asymmetry is then tranformed into baryon asymmetry by the standard model sphaleron 
processes violating (B + L)0. When the right-handed neutrinos are produced through scatterings (mainly inverse 
decays) in the thermal plasma this scenario is known as “thermal leptogenesis”. 

In the thermal leptogenesis mechanism, it is usually assumed that the lepton asymmetry is generated by the 
out-of-equilibrium decays of the lightest right-handed neutrino, N±. The asymmetry generated by the two heavier 
right-handed neutrinos, N 2 and iV 3 that we take as hierarchically heavier than iVi, is usually assumed to be either 
negligible or efficiently washed-out by the IVi inverse decays 0- In these conditions only iVi generates a baryon 
asymmetry. The final value of the baryon asymmetry normalized to the entropy density, 773, is then related to the 
CP asymmetry, gq, in the decays of the lightest right-handed neutrino as 
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where the factor 8/23 is the fraction of the B — L asymmetry converted into baryon asymmetry by sphalerons. 
i^Ni/s is the number density of right-handed neutrinos normalized to the entropy density that in equilibrium is 
approximately equal to 0.2 /g* with g * ~ 230 the number of propagating states in the supersymmetric plasma. Finally 
Xi is the efficiency factor describing the fraction of the CP asymmetry that survives the washout by inverse decays 
and scattering processes. Now, for hierarchical right-handed neutrinos, it is possible to relate £\ to the mass of the 
lightest right-handed neutrino, the observed neutrino masses, and an effective leptogenesis phase as 
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where the maximum value for this asymmetry corresponds to fully hierarchical neutrinos with m 3 = rn atm ~ 0.05 eV, 
~ 0 and a maximal phase 8^. 

The baryon asymmetry is precisely measured from the cosmic microwave background as ?7 b MB = 0.9 x 10“ 10 jij. 
Thus, using Eqs. 0 and 0 with xi = 1 we obtain a lower bound on the mass of N x 
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However this lower bound represents a serious problem for many flavour models, specially those incorporating some 
grand unification symmetry. In these models we usuallvexpect that neutrino Yukawas are closely related to the 
up-quark Yukawa matrices. As show in the appendix of p (see also Ref. i) and explained below, in these models 
it is possible to obtain the mass of the lightest right-handed neutrino as 
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with Wij the rotation from the basis of diagonal jzl masses to the basis where the (left-handed) neutrino Yukawa 
matrix Y u Yg is diagonal, m Uk the left-handed neutrino masses that we take as hierarchical and we have approximated 
yiV 2 — rn u . From this expression it is straightforward to see that we can only obtain M\ > 10 9 GeV if 
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This implies that, barring accidental cancellations (or some specific textures as in Jgj), having a large enough N± is 
only possible if we have simultaneously mi < 4x 10 -6 eV, W 2 2 < 5 x 10 -4 (for m so i — 0.008 eV) and W 2 3 < 8 x 10 -5 
(for m atm ~ 0.05 eV). These conditions are usually not satisfied in most flavour models where the typical values for 
the lightest right-handed neutrino mass are close to 10 6 -10 7 GeV with W 32 of the order of mixing angle in solar 
neutrinos p, 113, [ll| ■ Therefore, apparently, these models can not produce a sufficient baryon asymmetry through the 
thermal leptogenesis mechanism. 

It is evident that this represents a very serious problem for all these flavour models, that are forced to abandon 
thermal leptogenesis and use other mechanisms to produce the observed baryon asymmetry 0 In this paper we 
present a solution to this problem. We show that when the flavour dependence of the lepton asymmetries and the 
erasure processes are correctly taken into account, it is still possible to generate a large enough baryon asymmetry 
in the decays of the second right-handed neutrino N 2 that then survives the washout due to the inverse decays of 
Ni. The basic idea behind this statement is that the asymmetries generated by N 2 in the different flavour channels 
are washed out by different N- t Yukawa interactions. In fact, the decays of the second right-handed neutrino N 2 
create different asymmetries in the r, y and e channels. As we will show these asymmetries do not mix through the 
thermal interactions with the plasma before the Ni inverse decays become effective. Therefore, it is evident that the 
different flavour channels will only be erased by the corresponding N\ —> HLi interaction and not by the sum over all 
lepton flavours as usually done |3, a [lj, liii llD, |l(j ■ Taking into account these effects, we prove here that in models 
with Y v — Y u , the lepton asymmetry in the r and /r channels is only mildly erased by Aq decays. Some of these 
ideas were already present in the literature. For instance, the possibility of generating the asymmetry in N 2 decays 
was previously discussed by P. di Bari in the one flavour approximation 0 ; an early discussion of flavour effects in 
leptogenesis, mainly in N- t decays, can be found in 0 , and flavour effects in the context of low scale leptogenesis 
models were recently studied in 1181 . In this paper we show that combining these elements it is possible to obtain 
successful leptogenesis in realistic models of fermion masses with up-quark-neutrino Yukawa unification. Therefore, 
we reopen the possibility of “thermal leptogenesis” in these models. 

In the following we give a general proof of this statement. Although this mechanism is equally valid in supersym¬ 
metric and non-supersymmetric models we only discuss the non-supersymmetric example. The supersymmetric case 
is completely analogous [Til . 

The general expression for leptonic Yukawa couplings and Majorana masses in a seesaw Standard Model is, 
£yuk = L T Y e e r Hi + L T Y u v g H 2 — A AIi'r, with Li,em,VRi , (* = e,y,r) the three generations of 
leptons and Y e 1 Y v and Ad, 3 x 3 matrices. In the basis of diagonal Y e and Ad, the neutrino Yukawa matrix is 
Y v = ■ Dy u • Vr, where Dy v is a diagonal mass of neutrino Yukawa couplings, that in the class of models dis¬ 

cussed here is Dy u = diag(yi, y 2 , y 3 ) ~ diag(m u /v 2 ,m c /v 2: mt/v 2 ). The effective Majorana mass for the left-handed 
neutrinos is obtained through the seesaw mechanism as 

m v = v\Y v -M~ l -Y vT = U- D mv ■ U T , (6) 


where the mixing matrix U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix that is sufficiently known from 
neutrino experiments, and D mv = diag(?ni, m 2 , m 3 ) the three neutrino masses. At the moment, we know only two 
mass differences, Am 2 tm ~ 2.7 x 10” 3 eV 2 and Am^ tm ~ 7.0 x 10 -5 eV 2 . Although in the following we assume 
hierarchical spectrum with normal hierarchy, other cases can be easily treated in the same way. 

Now we define the hermitian matrix A = v 2 yf Ad^^Ad -1 6} and, in the basis of diagonal Y V ^Y V , from Eq. © 
we have 
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with A = Vl m u V^ = W D mi/ W T . From here the matrix A can be written as 
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Given the strong hierarchy in the neutrino-Yukawa eigenvalues (similar to the hierarchy in the up-quark masses) we 
have in practice always that |Ai | 2 | A 2 1 2 , | A 3 1 2 . For instance, to change this hierarchy in the simplified case where we 

take W 13 — 0 and j/ 1 / 7/3 — 0 would require that |An| < \A 22 \y\/y 2 . In these conditions it is very easy to understand 
that the largest eigenvalue of the matrix A, corresponding to v 2 yf/\Mi\ 2 , will be given in very good approximation 
by the (1,1) element of this matrix. Therefore the lightest right-handed Majorana neutrino mass is given by: 
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From this expression we can estimate the mass of the lightest right-handed neutrino. If the observed neutrino mixings 
are not present in the neutrino-Yukawa matrices and come mainly from the seesaw mechanism itself we have that 
Vl — 1 and W ~ Upmns- Then we have that \M\\ ~ 2 rn^/m so \ ~ 1 x I0 6 GeV. Only in the case W ~ 1 we can 
have \Mi\ ~ m^/mi = 10 s (4 x 10 _ 5 eV/TOi) GeV. Therefore, typically Ni is too light to generate a sizeable baryon 
asymmetry. In the same way the associated eigenvector to this eigenvalue is given by Ai ~ Ai/An. Using this 
eigenvector we can also obtain the Yukawa couplings of Ni in the basis of diagonal Y e and Ai as 
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where we used that (U L )*- A a = J^ k Uj k m Vk Wi k . 

At this point we can already discuss the behaviour of a lepton asymmetry created by the out-of-equilibrium decays 
of the next-to-lightest right-handed neutrino N 2 . The asymmetry in the different flavour channels (a = e,/x, r), 
assuming Mi <C M 2 -C M 3 , is given by ; 13J 
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Owing to the hierarchy in the neutrino Yukawa matrices and the right-handed Majorana neutrinos, the first term 
in Eq. CD usually dominates and therefore we take it to estimate the size of the asymmetry. Now we have that, 
l/u 2 Im [Y* 2 ( m v ) ak Y fc2 ] = 1/M 3 Irn [Y* 2 Y a3 Y£ 2 Y k3 ] + l/Mi Im [Y a 2 Y 0 iY fc 2 Yfci]. Therefore, the numerator in e 2 includes 
the M 3 contribution to the {nii/) ak element. To maximize the asymmetry we assume this element to be of order 
m „ 3 = vTiatm and define the corresponding CP-violating phase as 8 Vl then we have 
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where we use that the hierarchy in the Yukawa elements is such that Y 3 j > Y 2 j , Y\ 3 . It is clear that asymmetry in the 
r channel can be as large as e 2 = 3 x 1CV 6 (M2/10 10 GeV). In the muon and electron channels the asymmetry will be 
further suppressed by IY 22 I/IT 32 I and |Y 12 1/1YA 2 1 respectively, and thus we can expect a smaller asymmetry specially 
in the electron channel | 22 |. 

These asymmetries are partially converted, through sphaleron processes, to a A a = Ip — L a asymmetry and, in 
part, washed away by N 2 inverse decays. The final lepton asymmetry after the decay of N 2 will be given by 
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with ri]\r 2 /s the number density of N 2 normalized to the entropy density[2^. The efficiency factor, \ 2 i th e fraction 
of the produced asymmetry that survives after the end of N 2 decays and inverse decays and it is approximately given 
by 
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with m 2 = v 2 Y a2 Y* 2 /M 2 13 and m* is the equilibrium mass equal to ~ 1 x 10 3 eV in the SM 13- The final values 
of this efficiency factor are model dependent and therefore we keep this factor as a parameter in the following. 

At temperatures slightly below M 2 we have three different asymmetries in the A r , A M and A e channels. These 
asymmetries remain basically unchanged from T ~ M 2 to T ~ Mi and, in particular, the different flavours do not 
mix. In fact, at this temperature, both N 2 and N 3 right-handed neutrinos have already decayed. Therefore the only 
possible Yukawa interactions in the plasma are those with the right-handed charged leptons and with Aq. Taking 
into account that the Yukawa interactions are in thermal equilibrium at temperatures below ~ 10 14 GeV, it is 
more convenient to use the basis of diagonal charged lepton Yukawas to include these interactions in the thermal 
masses. Thus, only the Yukawa interactions with Ni can change flavour, although iVi are very rare in the plasma up 
to temperatures of the order of Mi when they reach the equilibrium density |l5|,|l6j|. Therefore, we have to discus the 
washout of the different asymmetries by the iVi decays and inverse decays. First, it is easy to see that this washout 
will be exponential and not linear if the iVi are thermally produced at temperatures close to Mi. From Ref. 13 and 
taking e\ — 0, the Boltzman equation for this case would be j3| : 


d?7A a 

dz 


-l z 2 e- z 

4~ 


. 7T TO 1 

1 + 77 “ - A a b VA b , 

2 TO* 


/ — 0.86 0.1 0.1 \ 

with A(T < 10 9 GeV) ~ 0.06 -0.65 0.017 , 

\ 0.06 0.017 -0.65 J 


( 15 ) 








4 


with m° = Wj ^ai^ai/^i takes into account that a given asymmetry A a can only be erased by the corresponding 
lepton flavour, z = M\/T and m* ~ 1 x 10 ~ 3 eV[25}. Now we can use Eqs. © and m and we see that 
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where m 3 = m a tm and m 3 = m so i and we have used bimaximal mixings in the PMNS matrix to make an estimate. 
The final values of m“ depend on the matrix W. In case the case the mixings in the neutrino Yukawa matrix are 
small, we have that W ~ U and so 
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Using these values in the case 77 ^“ ~ 77 ™ 77 ^*, Eq. IT5ll can be easily integrated (taking ?7a t {z) — r/A^ (z) , r]A e ( z ) ~ 0 

and A tt + A Tll = 0.63) and we obtain 
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If now we take the central value from the solar mass difference as m so 1 = 0.008 eV, we have that the muon and tau 
asymmetries are only erased by a factor exp(—1.5) — 0.22 and hence they can survive the washout from N 3 decays. 
This should be compared with the washout that we would obtain in the single flavour approximation, that would be 
~ e -o .75 E.™"/" 1 * = e -0 - 75 m noi/m, _ 2.4 x 10 -3 . Now the final asymmetry in the r channel would be 
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where we have used that t]a t — and rj a c — 0. Thus in this case, it could be possible to generate a sufficient 
baryon asymmetry from the p, and r channels even for efficiency factors x 2 an d X 2 order 0.15 and M 2 ~ 10 10 GeV. 
For smaller efficiency factors we could still have successful leptogenesis with a somewhat heavier M 2 . This situation 
is typically found for instance in models with non-Abelian flavour symmetries. In this paper we have only presented 
a simple estimate to show that this mechanism can generate a large enough baryon asymmetry. A full computation 
in a model with an 517(3) flavour symmetry and spontaneous CP violation is now in progress 0, 11, l21| . 

From Eq. ITHl) and taking into account that An = W 2 k m Vk we can see the washout by N\ can only be stronger 
than the result obtained in Eq. if the contribution from the atmospheric neutrino mass is sizable (assuming that 
TOi <C m so i). This requires basically that VE 13 TO 3 > fVi 2 rn 2 . Clearly this possibility can not be excluded, but taking 
into account that m 2 /m 3 ~ 1/6, W = Vl ■ U and that U\ 2 ~ l/x/2 and U 13 < 0.2 it requires large left-handed mixings 
in the neutrino Yukawa matrix. This could be possible in models based in Abelian flavour symmetries or discrete 
symmetries, although it is clear the final result is highly model dependent and the different models must be studied 
in detail. 

We have shown that taking correctly into account the lepton flavour dependence of the CP asymmetries and 
washout processes, it is possible to obtain successful thermal leptogenesis from the decays of the second right-handed 
neutrino. The asymmetries in the muon and tau-flavour channels are then not erased by the inverse decays of the 
lightest right-handed neutrino, N±. Therefore “thermal leptogenesis” is still viable in models with a strong hierarchy 
in the right-handed Majorana masses that is typically the case in models with up-quark-neutrino Yukawa unification. 
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